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Sliding parts in nanosystems such as Nano 
ElectroMechanical Systems (NEMS) and 
nanomotorsi - — , increasingly involve large speeds, 
and rotations as well as translations of the 
moving surfaces; yet, the physics of high speed 
nanoscale friction is so far unexplored. Here, by 
simulating the motion of drifting and of kicked 
Au clusters on graphite a workhorse system 
of experimental relevanceiSr-i^- we demonstrate 
and characterize a novel "ballistic" friction 
regime at high speed, separate from drift at 
low speed. The temperature dependence of the 
cluster slip distance and time, measuring friction, 
is opposite in these two regimes, consistent with 
theory. Crucial to both regimes is the interplay 
of rotations and translations, shown to be cor- 
related in slow drift but anticorrelated in fast 
sliding. Despite these differences, we find the 
velocity dependence of ballistic friction to be, 
like drift, viscous. 

The friction of adsorbed molecules and mobile clus- 
ters thermally diffusing on surfaces is gradually moving 
from a theoretical concept to reality as shown for exam- 
ple by recent Helium-3 spin-echo experiments^. As well 
known in Brownian motion^, everything that diffuses 
will slowly drift and exhibit viscous friction under a weak 
force; we are however much more ignorant about the fric- 
tional behaviour to expect when the adsorbatc is forced 
to move at a higher externally imposed speed. High 
speed friction has long been known in the traditional 
macroscopic context^, and also in AFM experiments 17 
on asperity dominated rough surface friction, with highly 
case-specific outcomes. Here we address the simpler and 
more fundamental question of speed dependent friction 
during sliding motion of nanosized objects on atomically 
flat surfaces, a situation more likely to yield a result of 
generic validity. To gain experience, we undertake to ex- 
plore the problem by simulations. Adsorbed Au clusters 
on graphite^ are known to be thermally mobile even at 
room temperature^, and make an ideal test case. Ex- 
perimentally, the low speed Au cluster friction is in prin- 
ciple accessible in Quartz Crystal Microbalance (QCM) 
experiments^—, while high speed friction could be re- 
alized in AFM/STM tip based setups, where it is possi- 
ble to pursue kicking techniques, either mechanica l 22 i 23 
or electrical. Anticipating experiment, we simulate by 
realistic molecular dynamics (MD), the speed-dependent 



friction of small Au clusters (typically 250-500 Au atoms, 
as in Fig. |1} ona fully mobile graphite surface (see Meth- 
ods and Supplementary Information). 

Canonical MD is initially used to simulate the force- 
free cluster diffusion as a function of temperature T. 

The Au cluster center-of-mass (CM) coordinate R(t) 
is found to diffuse positionally on graphite, with re- 
sults that are close both to experiment and to previous 
simulations^— . In addition we observe that the two di- 
mensional CM positional random walk is closely linked 
to another conjugate, and so far unexplored, one dimen- 
sional random walk executed by the cluster orientation 
angle 6(t) (Fig. HJa)). For relatively long time lapses 
there is good angular registry between the cluster's (111) 
contact facet and the graphite substrate (see bottom 
panel (I) in Fig. [IJ, and both translations and rotations 
are stuck- only small oscillations occur. Once a ther- 
mal fluctuation depins the registry angle 9(t), the clus- 
ter rotates out of commcnsurability with the substrate, 
(bottom panel (II) of Fig. QJ. Now free to slide, the ro- 
tated cluster glides long-distance over the surface, with 
large jumps which dominate diffusion (see Supplemen- 




FIG. 1: Simulated gold cluster on graphite substrate. 

Top panel: A truncated octahedron AU459 cluster sliding with 
its 36-atom (111) facet (area 177 A 2 ) over a mobile graphite 
substrate. Bottom panel: In-registry (I) and out-of-registry 
(II) geometries for the cluster's (111) contact facet and the 
graphite substrate. 
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FIG. 2: Thermal diffusion of deposited Au clusters on 
graphite. Panel (a): Force free positional and angular ther- 
mal evolution at T = 700 K. Pinned states corresponds to par- 
allel orientations of the (111) cluster facet relative to the sub- 
strate (see panel (I) in Fig. [TJ. Note the correlation between 
positional and angular stops and jumps. Inset: T-dependence 
of diffusion coefficients (dots), obtained from the first deriva- 
tive of the mean-square-displacement for long microcanonical 
simulations (> 100ns). Both positional and angular diffu- 
sion fit a thermally activated form D = Do exp (— Eo/ksT) 
(continuous lines). Panel (b): Predicted temperature depen- 
dent QCM diffusional slip time of a drifting AU250 cluster on 
graphite, as expected by Einstein's relation r = mDa/kBT 
with D(T) = Doexp(—Ed/kBT), (experimental diffusivity 
values used**, D = 10 3 cm 2 s _1 , E d = 0.5 eV. 



tary Information), as recognized long ago by Luedtke et 
al. and by others^— . We find here that both positional 
and angular random walks are well fit by R 2 ~ AD at 
and 9 2 ~ IDgt. For variable temperature, both dif- 
fusion coefficients are thermally activated (see inset in 
Fig. Ha)), D = D exp(-E /k B T) with D 0R = 0.02 
cm 2 /s, Dqs = 0.77 x 10 14 rad 2 /s with activation barriers 
Eor = 0.351 eV and E e = 0.509 cV respectively, un- 
derlining the controlling role of angular depinning over 
positional diffusions*. In this respect, it would be inter- 
esting if like positional diffusion in He3 spin-echo 1 -^, the 
angular diffusion of molecules too could be accessed by 
some experimental technique. 

Under a weak applied sliding force, such as realized, 
e.g., in QCM experiment o 19 ' 20 , the freely diffusing clus- 



ters or molecules will respond by drifting. Initial simu- 
lations in presence of a weak applied force (see Supple- 
mentary Information), show that during drift the strict 
correlation seen in diffusion between angular depinning 
and positional jumps is preserved. In QCM one mea- 
sures the slip time r = m/i, with typical detectable values 
in the order of 10 _9 -10 -8 sec. Full fledge drift simula- 
tions would of course require prohibitively long times; 
luckily here we expect drift and diffusion to obey linear 
response, connecting the positional drift mobility fj, to 
the viscous diffusion coefficient Dr through Einstein's 
relation Dr = /iksT. Together, linear response and 
activated cluster diffusion predict a frictional slip time 

r = ^gfj^exp f-fkf). Using that, data of Bardotti 

et aZ3 with a room temperature diffusion coefficient 
D ~ 10 -5 cm 2 /sec for AU250 nanoclusters on graphite, 
translate to a slip time r ~ 10 -11 sec, probably too 
short for experimental detection, but rapidly growing to 
measurable values ~ 10 -9 or larger at 400 K and above 
(Fig. HJb))2i. Actual QCM experiments would be desir- 
able to check this diffusional friction of Au clusters on 
graphite. 

We come now to high speed nanofriction, our main 
task. Here, the linear response theory connection to dif- 
fusion no longer holds and simulations are essential. A 
cluster receiving a large side kick will start off with a large 
speed vq and slide inertially, its center of mass slowing 
down with some as yet unknown time dependence ve- 
locity v(t). We will refer, for lack of a better name, to 
this frictional regime as "ballistic" , in analogy with the 
damped motion of a projectile. At high speed, the static 
substrate corrugation is averaged out, and ballistic fric- 
tion will arise from collisions of the fast sliding cluster 
against the dynamic excitations of the medium. The de- 
fect free graphite surface has instantaneous and random 
thermal corrugations, whose magnitude < it 2 >oc k B T at 
high temperature. The number of collisions with these 
thermal excitations experienced by the fast cluster per 
unit time should be proportional to speed; and so there- 
fore should the energy damping in linear response. Sup- 
posing therefore provisionally a viscous slowdown (same 
as for diffusive drift) dv/dt = —jv + 77, (here 7 is a fric- 
tion coefficient and 77 a random force), friction would 
imply the exponential decay of cluster speed with time. 
Thermal corrugations should provide a frictional damp- 
ing 7b = Ffriction/v oc k B T, and thus a "ballistic slip 
time" t;, = ^ and a ballistic slip distance A&, (the latter 
defined as the distance reached by the kicked cluster be- 
fore major slowdown, equal to Ah = «T(,(1 — e -1 )) oc j^pf- 
We therefore expect that the ballistic slip time and dis- 
tance will drop with temperature, contrary to diffusive 
sliding, where instead they grew (Fig. Efb)). We note 
however that the significance of this slip time is quite 
different from the diffusive one. 

We can directly compare this speculative scenario with 
actual simulations of high speed cluster sliding. The MD 
ballistic simulations were initialized by first thcrmalizing 
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the cluster-surface system at zero speed and temperature 
T (for about 100 ps), then removing temperature control 
while giving the cluster a large initial sliding speed 
(with a typical kickoff value vq = 100m/s, see Methods 
and movies in Supplementary Information). 

Figure [3] portrays the (X, Y) cluster CM trajectories 
of kicked AU459 for increasing temperatures. The cluster 
motion is initially ballistic and straight with a slowdown 
reasonably fit (despite a sizable error caused by irregu- 
lar motion) by the viscous exponential law (Fig. [H top 
panel). The straight trajectory tract in Fig. [3] provides a 
visual measure of the slip distance A;,, with a correspond- 
ing slip time which drops with T, close to the expected 
\/T behaviour (inset of Fig. |3J) ■ The straight ballis- 
tic tract generally ends with a sharp turn, followed by 
a crossover to a curved, slower and increasingly erratic 
trajectory. 

The role of rotations is important in ballistic fric- 
tional slowdown; in fact cluster rotations are closely con- 
nected with translations, as was the case in diffusion 
and drift. When no angular momentum is initially im- 
parted, the cluster angular dynamics is relatively mod- 
est in the straight initial part of trajectories. However, 
as soon as the cluster-substrate orientations match, the 
CM translational energy exhibits a sudden drop in fa- 
vor of an increase of the cluster angular speed, subse- 
quently reversing back to translational, and so on. To 
make this "push-pull" anticorrelation between transla- 
tions and rotations more visible we simulated clusters 
with larger moments of inertia and contact facets. MD 
data for AU2406 in bottom panel of Fig. H] show a CM 
slowdown consisting of tilted plateaus for out-of-registry 
cluster-graphite configurations and sharp drops at acci- 
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FIG. 3: (X,Y) ballistic trajectories of a kicked cluster 
at different temperatures. Top view of center of mass tra- 
jectories of a Au45g cluster with initial velocity v = 100 m/s 
at different temperatures. Note the initial straight ballistic 
tract, followed by a broadly curved crossover region. The 
inset shows the ballistic slip time rt, obtained from an expo- 
nential fit of the length of the straight slowdown trajectory 
portion. The decrease of rt, with temperature (oc j 1-0 - 82 ) j s 
close to the expected 1/T. 
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FIG. 4: Slowdown dynamics of kicked Au clusters. Top: 
Slowdown and crossover between ballistic to diffusive friction 
of a small cluster. CM speed, displacement, and cluster rota- 
tional angle are shown. An exponential velocity slowdown is 
superposed to the velocity data. At t = lOOps, the thermal- 
ized cluster is kicked with a high speed vo — lOOm/s. Bottom: 
detail of ballistic slowdown of a larger AU2406 cluster, high- 
lighting clearer sharp drops at exact matching angles of the 
cluster contact (111) facet with the graphite substrate. 



dentally perfect matching angles. Connections between 
oricntational mismatch and easy sliding was underscored 
in previous force free diffusion simulations^. In ballis- 
tic slowdown the speed plateaus degrade gently, corre- 
sponding to "superlubric" incommensurate friction, well 
known, e.g., for graphite flakes on graphite^; here, this 
superlubric dynamical state is favored by the fast sliding 
motion, preventing significant angular pinning. As the 
trajectories cross over from straight to curved, the clus- 
ter angular and positional evolutions become concerted 
with some analogy with graphite flake behaviour—. The 
turning point from a straight to curved trajectory marks 
the beginning of a broad crossover between initial bal- 
listic motion and final diffusive drift (see Supplementary 
Information). At the end of the crossover region (whose 
extension depends on inertia and on T) the, now corre- 
lated, stop-and-go motion typical of diffusion makes its 
appearance, the cluster-substrate angular registry again 
corresponding to stops. This strong translation-rotation 
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interplay is reminiscent of the macroscopic theory on slid- 
ing disk o 26 ' , where the friction force and torque are in- 
herently coupled, although differently from here. To shed 
further light on these aspects, we compare the dynamics 
of clusters initialized with a simple linear velocity Vo to 
that of clusters initialized simultaneously with a linear 
speed vq and an angular velocity wo (see movie in Sup- 
plementary Information) . Their kinetic energy evolutions 
are respectively shown in top and bottom panels of Fig. 
[5j On average, both v(t) and u(t) decay in an approx- 
imately exponential manner. The platcau-drop-plateau 
CM motion pattern here disappears, as rotation prevents 
the touching surfaces from locking together. Fluctuations 
during the ballistic slowdown show a clear anticorrelation 
between angular and translational speeds, as highlighted 
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FIG. 5: Time decay of kinetic energies of purely sliding 
and of sliding/rotating ballistic Au clusters. Ballistic 
time dependence of translational (red), rotational (blue), and 
total (green) kinetic energies of the AU459 cluster at 300 K, ini- 
tialized with a linear CM velocity vo = 100m/s(l A/ps) and 
zero initial angular momentum [top], or with an additional 
large initial rotational speed 1^0 =0.1 rad/ps [bottom]. Note 
the positional-orientational anticorrelation, as CM slowdowns 
trigger angular speedups. In the rotating case, time decays of 
linear, angular and total kinetic energies are fit by exponen- 
tials. The inset highlights the anticorrelation between cluster 
rotation and translation, obtained by subtracting the smooth 
exponentials from v(t), and u)(t). 



in the figure inset. This is just the opposite of the correla- 
tion found in diffusion and in drift motion, where angular 
and positional stops and goes coincide. By analysing the 
separate translational and rotational kinetic energies we 
find that in fast rotating clusters the anticorrelation be- 
tween translational and rotational speed drops is directly 
connected with overall kinetic energy near conservation 
on a short time scale, resulting in an energy push-pull 
between the two degrees of freedom. As the skidding 
cluster collides with a surface thermal bump it suddenly 
slows down. The bump is typically not central, and the 
hit causes a rotation speedup, with a temporary push of 
kinetic energy from translation to rotation. The higher 
rotation speed now increases the chance to symmetri- 
cally hit another bump, whereby kinetic energy is pulled 
back to translation; and so on. The total translation 
plus rotation kinetic energy decreases as a result, as the 
bottom panel of FigJS] shows, much more smoothly and 
exponentially than either of the two components, much 
bumpier. Ballistic translation-rotation anticorrelation is 
also evident in the purely translational data of FigH] or 
FigEJ top panel. In general, we note that the total kinetic 
energy decay is well fit by a simple exponential, imply- 
ing that the overall energy damping is indeed reasonably 
described by dE/dt oc — E, confirming that the ballis- 
tic friction of high speed translating-rotating clusters is 
again viscous, as in low speed diffusive friction. The rea- 
son however seems quite different: simply, in ballistic 
friction the number of thermal asperities encountered by 
the slider must to be proportional to speed. 

In summary, the frictional sliding of diffusive clusters 
(or large molecules, whose behavior should be similar) 
illustrates two distinct nanofriction regimes - a standard 
diffusional one at low speed, and a new ballistic one at 
high speed - where the effect of temperature on friction 
is opposite. Angular motion and rotations cannot be 
ignored and are crucial to both diffusive and ballistic 
regimes. The interplay between translations and rota- 
tions, besides turning ballistic trajectories from straight 
to curved, controls the crossover between the two sliding 
regimes. Translational and rotational stops and goes, 
strictly correlated in diffusion and drift, become anticor- 
related in fast ballistic sliding. Ballistic friction, although 
quite different from diffusive friction and naturally erratic 
in small clusters, appears still viscous in its speed depen- 
dence, a useful result whose validity for this and other 
systems was not discounted in principle, and which it 
would be interesting to pursue and test experimentally. 
Some evidence, we note, may already be present in fast 
AFM friction dataii. Angular diffusion might conceiv- 
ably be addressed by labeling a cluster periphery, and 
the diffusive drift of clusters should be directly accessible 
in QCM experiments 2 ^.. Extending the well developed 
tip-based studie a 22 ' 23 , the possibility to kick clusters bal- 
listically either mechanically or electrically through e.g., 
a charge pulse, seems entirely open to experimental veri- 
fication. Finally, advances in fundamental understanding 
of ballistic nanofriction should impact practical applica- 
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tions where the relative speeds are not small. 

METHODS Au fee clusters were generated and their 
shape optimized with an empirical EAM potential^. 
Among several, a truncated octahedron A11459 cluster was 
selected for most sliding simulations (Fig. [1]). Dynamics 
of a two- layer and a three- layer fully mobile graphite sub- 
strate was considered and described with Tersoff intra- 
layer potential^ and Kolmogorov-Crespi RDPO inter- 
layer interactions^. Interatomic Au-C interactions were 
parametrized by a standard Lennard- Jones form, with 
parameters yielding realistic corrugation barriers of the 
graphite substrate against cluster sliding^ (see Supple- 
mentary Information). 

Temperature equilibration in force free MD simulations 
was enacted by a standard Langevin thermostat. All sub- 
sequent frictional simulations were carried out after re- 



moving the thermostat, and carefully checking that the 
frictional energy released into the substrate was small 
enough to cause a negligible substrate temperature in- 
crease. Unlike previous force free simulations^, where a 
rigid substrate sufficed to study diffusion, the fully mo- 
bile graphite substrate employed here is essential to study 
sliding friction. We checked that two or three mobile 
graphite layers gave sufficiently similar results to mimic 
the infinite substrate. 

Finally, in kicked simulations the initial cluster velocity 
was given a large value of 100 m/s where the physics of 
the ballistic regime is extractable for a substantially lower 
computational cost than for lower speeds. Results show 
that crossover from ballistic to diffusive occur between 
10 and 1 m/s. 
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